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Abstract 

We prove that the electromagnetic fields in dielectric media whose 

susceptibility follows a fractional power-law dependence in a wide fre- 

CN ■ quency range can be described by differential equations with time 

CN ! derivatives of noninteger order. We obtain fractional integro-differential 

ON . 

OO ■ equations for electromagnetic waves in a dielectric. The electromag- 

^^ , ation. The fractional integro-differential equations for electromagnetic 



Ph. 



% 



waves are common to a wide class of dielectric media regardless of the 
type of physical structure, the chemical composition, or the nature of 
the polarizing species (dipoles, electrons, or ions). 
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1 Introduction 

Debye formulated his theory of dipole relaxation in dielectrics in 1912 fT]. A 
large number of dielectric relaxation measurements show that the classical 
Debye behavior is very rarely observed experimentally [21 El S]. Dielectric 
measurements by Jonscher for a wide class of various substances confirm that 
different dielectric spectra are described by power laws [21 [3] . 

For the majority of materials, the dielectric susceptibility in a wide fre- 
quency range follows a fractional power-law called the universal response 
[21 [3]. This law is found both in dipolar media beyond their loss peak fre- 
quency and in media where the polarization arises from movements of either 
ionic or electronic hopping charge carriers. It was shown in [S] that the 
frequency dependence of the dielectric susceptibility x{uj) = x'{^) ~ ^x"{.'^) 
follows a common universal pattern for virtually all kinds of media over many 
decades of frequency, 

XV)~^"'\ x'V)~u;"-\ (u;»u;p), (1) 

and 

x'(0)-xV)~^™ x"M~c^^ (c^«cup), (2) 

where x'(0) is the static polarization, < n, m < 1, and Up is the loss peak 
frequency. We note that the ratio of the imaginary to the real component 
of the susceptibility is independent of frequency. The frequency dependence 
given by equation ([T]) implies that the imaginary and real components of the 
complex susceptibility at high frequencies satisfy the relation 

|ii = cot(-). (.».,). (3) 

Experimental behavior ([2|) leads to a similar frequency-independent rule for 



the low- frequency polarization decrement, 

The laws of universal response for dielectric media [21 [3] can be described 
using fractional calculus [6] . The theory of integrals and derivatives of non- 
integer order goes back to Leibniz, Liouville, Riemann, Grunwald, and Let- 
nikov [6]. Fractional analysis has found many applications in recent studies 
in mechanics and physics. The interest in fractional equations has been grow- 
ing continuously during the last few years because of numerous applications. 
In a short time, the list of applications has becomes long (see, e.g., [3 [HI [9]). 
In Refs. [ini [m [121 [13] , fractional calculus has used to explain the nature of 
nonexponential relaxation, and equations were obtained containing operators 
of fractional integration and differentiation. 

Here, we prove that a fractional power-law frequency dependence in a time 
domain gives integro-differential equations with derivatives and integrals of 
noninteger order. We obtain fractional equations that describe electromag- 
netic waves for a wide class of dielectric media. The power laws of Jonscher 
are represented by fractional integro-differential equations. The electromag- 
netic fields in the dielectric media demonstrate universal fractional damping. 
The suggested fractional equations are common (universal) to a wide class 
of materials regardless of the type of physical structure, the chemical com- 
position, or the nature of the polarizing species. 



2 Fractional equations for universal laws 

We consider Eqs. ([1]) and ([3]). For the region u ^ Wp, universal fractional 
power law ([1]) can be presented in the form 

xM = XaM"", (0<a<l) (5) 

with some positive constant Xa and a = 1 — n. Here, 

(icj)" = |a;|"exp{ia7rs(?ri(w)/2}. 

It is easy to see that relation ([3]) is satisfied for ([5]). 
The polarization density P(t,r) can be written as 

P(t, r) = 7-' {P{uj, r)) = e,F-' (x(c^)E(u;, r)) , (6) 

where P(a;,r) is the Fourier transform T of P(t, r). Substitution of ([5]) into 
dS]) gives 

We note that the Fourier transform of the fractional Liouville integral 
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r(a) i_„ (t - t'Y-' 

is given by the relation (see Theorem 7.1 in [6] and Theorem 2.15 in [14j): 

where < Re{a) < 1 and f(t) e Li{R), or 1 < p < 1/Re{a) and f(t) € 
Lp{R). 

Using the fractional Liouville integral and the fractional power law ([5]) 
for x(a;) in the frequency domain, we obtain 

P(t,r) = £oXa(/+E)(t,r), (0 < « < 1). (7) 



This equation shows that the polarization density P(t, r) in the high-frequency 
region is proportional to the fractional Liouville integral of the electric field 
E(t,r). 

We consider Eqs. ([2]) and (j4]). For the region uo <^ Up, universal fractional 
power law (|2]) can be presented as 

X{u) = x(0) - x^i^ujf, (0 < /3 < 1) (8) 

with some positive constants X/3; x(0), and /3 = m. It is easy to prove that 
equation (jl]) is satisfied. 

We note that the Fourier transforms of the fractional Liouville derivative 
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where k — 1 < /3 < k, are given by the formula (see Theorem 7.1 in [6] and 
Theorem 2.15 in [H]): 

where < i?e(/3) < 1 and /(t) G Li(R), or 1 < p < l/Re{P) and f{t) e 



Using the definition of the fractional Liouville derivative and fractional 
power laws ([8]), we can represent polarization density ([6]) in the form 

P(t, r) = £oX(0)E(t, r) - £oX^(D^E)(t, r), (0 < /3 < 1). (9) 

This equation shows that the polarization density P(t, r) in the low- frequency 
region is determined by the fractional Liouville derivative of the electric field 
E(t,r). 

Relations ([7]) and ([9]) can be considered universal laws. These equa- 
tions with integro-differentiation of noninteger order allow obtaining frac- 
tional wave equations for electric and magnetic fields. 



3 Universal electromagnetic wave equation 

Here, we obtain fractional equations for electromagnetic fields in dielectric 
media. Using the Maxwell equations, we obtain 

''m^ + -^^ + - {grad d^vE - V^E) + -^ = 0. (10) 

For u ^ Up, the polarization density P(t,r) is related with E(t, r) by 
equation ([7]). Substituting ([7]) in f lTO|) . we obtain the fractional differential 
equation for the electric field 

i^^ + ^(Dl-m, r) + {gradd^vE - V^E) = -/.^, (11) 

where < a < 1, and v"^ = l/(eo/^)- We note that divE ^ for p(t,r) = 0. 

For u ^ Wp, the fields P(t, r) and E(t,r) are related by Eq. ([9]). In this 
case, equation ( TTOl) becomes 



t//D Co (J o CJT/ 



where 

,2 _ 1 . _ X/3 



■y^ 



£o/u[l + x(0)]' '^ 1 + X(0)' 
Equations fITT]) and flT^ describe the time evolution of the electric field in 

dielectric media. These equations are fractional differential equations [H] 

with derivatives of the orders 2 — a and 2 + /3. 

Using the Maxwell equations, we obtain the equation for the magnetic 

field 

^^ = V^B{t, r) + -—curlP{t, r) + -curlUt, r). (13) 

In experiments, the field B(t, r) can be presented as B(t, r) = for t < 0, 
and B(t, r) ^ for t > 0. For u ^ Up, the polarization density P(t, r) is 
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related to E(t, r) by equation ([7]), which leads to the fractional differential 
equation for magnetic field in the form 

J_9^Bg^ + ^ (oA'-"B) (t,r) - V^B(t,r) = /.c«Hj(t,r), (14) 

where < a < 1, f^ = l/(£:oAi), and o-D^~" is the Riemann-Liouville deriva- 
tive [13] on [0, oo) such that 

For w <^ cjp, we obtain 
J_9^Bg^ - ^ (oA'^'b) (t,r) - V'B{t,r) = /icnr/j(t, r), (15) 

where < /9 < 1, and 



,2 



1 X/3 



^/3 - .„„ri ^.-.rmr ^/^ 



£o/i[l + x(0)]' '^ l + x(0)' 
Equations (iT4l) and ( IT5l) are fractional differential equations that describe the 
magnetic field in dielectric media and demonstrate a power-law relaxation. 
They can be written in a general form. Such a general fractional differential 
equation for the magnetic field has the form 

(oA°B)(t,r) - Ai (oA^b) (t,r) - X2V^B{t,r) = f(t,r), (16) 

where 1 < /3 < a < 3. The curl of the current density of free charges 
is regarded as an external source: f(t, r) = (1X2 curl i{t,r). Equation (fT6!) 
yields Eq. ([14]) for a = 2, 1 < /3 < 2, and 



Ai = -Xa, X2 = v'^ = l/(eoyu). 
Equation (TT5|) can be written in form (TT6!) for2<a<3,/3 = 2, and 

1 _ l + x(0) , _ ^} _ -1 
M — — — 5 ^2 — 



a/3 X/3 «/3 ^OAiX/3 



An exact solution of Eq. f ll6p can be written in terms of Wright functions 
using Theorem 5.5 in [2]. We note that Wright functions can be represented 
as derivatives of the Mittag-Leffler function Ea,j3[z] (see p3]). Solutions of 
equation ( [T6|l describe the fractional power-law damping of the magnetic 
field in dielectric media. An important property of the evolution described 
by the fractional differential equations is that the solutions have fractional 
power-law tails. 

4 Conclusion 

We have prove that the electromagnetic fields and waves in a wide class of 
dielectric media must be described by fractional differential equations with 
derivatives of the order 2 — a and 2 + /3, where < a < 1 and < /3 < 1. The 
parameters a = 1 — n and (5 = m are defined by the exponents n and m in the 
experimentally measured frequency dependence of the dielectric susceptibil- 
ity, called the universal response laws. An important property of the dynam- 
ics described by fractional differential equations for electromagnetic fields is 
that the solutions have fractional power-law tails. The suggested fractional 
integro-differential equations for the universal electromagnetic waves in di- 
electrics are common (universal) to a wide class of media regardless of the 
type of physical structure, the chemical composition, or the nature of the 
polarizing species (dipoles, electrons, or ions). 

We note that the differential equations with derivatives of noninteger 
order proposed for describing the electromagnetic field in dielectric media 
can be solved numerically. For example, the Grunwald-Letnikov discretiza- 
tion scheme [6] is used for numerically model the electromagnetic field in 
dielectrics described by fractional differential equations. For small fraction- 



8 



ality of a (or /3), an e-expansion [15] in the small parameter e = a (or 
e = 1 — /3) can be used. We note that a possible physical interpretation of 
fractional integrals and derivatives can be connected with memory effects or 
fractal properties of media (see, e.g., [Ml [T7]). 
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